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Abstract 



In this article, we introduce a general theoretical framework to analyze non-consistent 
approximations of the discrete eigenmodes of a self-adjoint operator. We focus in particular 
^Nj ' on the discrete eigenvalues laying in spectral gaps. We first provide a priori error estimates on 

the eigenvalues and eigenvectors in the absence of spectral pollution. We then show that the 
supercell method for perturbed periodic Schrodinger operators falls into the scope of our study. 
We prove that this method is spectral pollution free, and we derive optimal convergence rates 
for the planewave discretization method, taking numerical integration errors into account. 
^ (— I Some numerical illustrations are provided. 
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1 Introduction 

This article is concerned with the numerical analysis of the computation of the discrete eigenmodes 
of a self-adjoint operator A, on an infinite dimensional separable Hilbert space H. The focus is 
particularly set on the eigenmodes corresponding to discrete eigenvalues located in spectral gaps. 
The main application we have in mind is concerned with perturbed periodic Schrodinger oper- 
£^ . ators of the form 

A := -A + V pcl + W, 
IT} , 

where A is the Laplace operator on L 2 (M. d ), V per a periodic function of L^ oc (R d ) with p = 2 if d < 3, 
p > 2 for d = 4 and p = d/2 for d > 5, and W £ L°°(M. d ) a perturbation of the potential going 
to zero at infinity. The operator A is self-adjoint and bounded from below on H := L 2 (R d ) with 
domain H 2 (R d ). Perturbed periodic Schrodinger operators are encountered in electronic structure 
theory, and in the study of photonic crystals. In the case of a perfectly periodic crystal (W = 0), 
the spectrum of the operator A := —A + V pcr is purely absolutely continuous, and composed 
of a union of intervals of R. It follows from Weyl's theorem [23] that the essential spectra of A 
and A are identical. On the other hand, when W ^ 0, some discrete eigenvalues may appear in 
the band gaps of the spectrum of A. The corresponding eigenmodes, which can be interpreted as 
bound states trapped by local defects, are difficult to compute for numerical methods can produce 
spectral pollution. 

In a general theoretical framework, the eigenvalues of A and the associated eigenvectors can be 
obtained by solving the variational problem 

find (V>, A) e Q(A) x R such that , . 

€ Q(4), a(V>,0 = AmftM), U 

where m(-, •) is the scalar product of %, Q{A) the form domain of A, and a(-, •) the sesquilinear 
form associated with A (see for instance [9]). 
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A sequence (X n ) n< =fj of finite dimensional approximation subspaces of Q{A) being given, we 
consider for all neN, the self-adjoint operator A\x n ■ X n — > X„ defined by 

V(^„) <Pn) € X n X X n , m (A\ Xn 1p n , <t>n) = a(lp n , <t>n)- 

The standard Galerkin method consists in approximating the discrete eigenvalues of the operator 
A by the eigenvalues of the discretized operators ^4|x„, the latter being obtained by solving the 
variational problem 

J find (ip n , A„) <E X n x R such that 

\ V0„ € X n , a(lp n ,4>n) = Km(lp n ,(f> n ). 

According to the Rayleigh-Ritz theorem [23], under the natural assumption that the sequence 
(X n ) ne m satisfies 

V^eQ(A), inf U- <t> n \\ Q{A) — ► 0, 

this method allows to compute the eigenmodes of A associated with the discrete eigenvalues located 
below the bottom of the essential spectrum. Indeed, if A is bounded below and possesses exactly M 
discrete eigenvalues Ai < A2 < • • • < Xm (taking multiplicities into account) lower than min<7 ess (A), 
where a css (A) denotes the essential spectrum of A, and if |A"| ^h;™ y are * ne eigenvalues of 
A\x nl it is well-known that 

vi<i<M, a; 1 1 A , . 

n— >oo 

The situation is much more delicate when one tries to approximate eigenvalues which are located 
in spectral gaps of A since 

VM < j < dim X n , A™ J, minCT css (A). 

When dealing with the approximation of discrete eigenvalues of A located in spectral gaps, the 
standard Galerkin method may give rise to spectral pollution: some sequences (A„)„ e pj, where for 
each n, X n G (J {A\x n ), may converge to real numbers which do not belong to the spectrum of 
A. Spectral pollution occurs in a broad variety of physical settings, including elasticity theory, 
electromagnetism, hydrodynamics and quantum physics [TJ [3J [TUl [5J , and has been extensively 
studied in the framework of the standard Galerkin method [1I1[5J[TT1[T1[T1[T1[I71[IH]. We refer 
to [8j [19] for an analysis of spectral pollution for perturbed periodic Schrodinger operators. 

On the other hand, few results have been published on the numerical computation of eigenmodes 
in spectral gaps by means of non-consistent methods, based on generalized eigenvalue problems of 
the form 

J find (i/j n , X n ) € X n x M such that 

\ V0„ e X n , a n (ip n , <p n ) = X n m n (ip n , <p n ), 

where for all n € N, a n (-, •) and m n {-, •) are symmetric bilinear forms on Q(A), a priori different 
from a(-, •) and to(-, •). 

In this article, we consider a general theoretical framework to analyze non-consistent methods 
for the computation of the discrete eigenmodes of a self-adjoint operator. After introducing some 
notation and definitions in Sections l2.1l and [2~2l we state our main result (Theorem l3.1[) in Section[3] 
Theorem 13.11 provides a priori error estimates on the eigenvalues and eigenvectors in the absence 
of spectral pollution. Its proof is given in Section [4] 

In Section [S] we show that the supercell method for perturbed periodic Schrodinger operators 
falls into the scope of Theorem 13.11 We prove that this method is spectral pollution free, and 
we derive optimal convergence rates for the planewave discretization method, taking numerical 
integration errors into account. The corresponding proofs are detailed in Section [6] and some 
numerical illustrations are provided in Section [7] 

2 Approximations of a self-adjoint operator 
2.1 Some notation 

Throughout this paper, % denotes a separable Hilbert space, endowed with the scalar product 
77i(-, •) and associated norm || • ||%, and A a self-adjoint operator on % with dense domain D(A). 
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We denote by Q(A) := Dd^j 1 / 2 ) the form domain of A and by a(, •) the symmetric bilinear form 
on Q{A) associated with A. Recall that the vector space Q(A), endowed with the scalar product 
(•. -)q(A), defined as 

W>, g Q(A), <f>) Q[A) := m(V, 4>) + m (jA] 1 / 2 ^, \A\ l / 2 <^ , 

is a Hilbert space; the associated norm is denoted by || • ||q(a)- 

Example 2.1. Perturbed periodic Schrddinger operators A := — A + V pcT + W are self-adjoint 
semibounded operators on H :— L 2 (Bi d ), with domain D(A) := H 2 (Bi d ) and form domain Q(A) := 
H\R d ). 

For any finite dimensional vector subspace X of H such that X C Q(A), we introduce the 
following notation 

• ix '• X '-t H is the canonical embedding of X into H; 

• i* x ■ W — ► X is the adjoint of ix, that is the orthogonal projection from H onto X associated 
with the scalar product m (•, ■); 

• A|x : -X" — > X is the self-adjoint operator on X defined by 

VW^)axI, m(A|xlM) = afoM); 

• IT^ : W W and II^' 4 '' : Q(A) — > are the orthogonal projections onto X for 
(H,m(-, •)) and (•, '>q(A))) respectively. 

We set 

■ :=<T{Af\a d {A), 

where <7d(^4) is the discrete spectrum of A, and where <r(A) is the closure of c(A), the spectrum of 
A, in K := EU{±oo}. A spectral gap of A is an interval (£~,E+) such that E~,£+ £(7 css (4)nl 
and (£-,£+) n ct css (A) = (which implies that Tr(l ( _ 00jE - ] (A)) = Tr(l [s+j0o) (A)) = oo). As 
usual, 1 b denotes the characteristic function of the Borel set Bet. The discrete eigenvalues of 
the operator A in a spectral gap (£~~, £ + ), if any, are isolated and of finite multiplicities, but can 
accumulate at E~ and/or E + [23] . 

Let us finally recall the notions of limit superior and limit inferior of a sequence of sets of 
complex numbers (see for instance [9]). 

Definition 2.1. Let (_E„)„ s n be a sequence of subsets o/C. 

• The set lim E n (limit superior) is the set of all complex numbers A G C such that there exist 

a subsequence (£ , „ f .)/ c eN of (E n ) n< zjq and a sequence (X nk )kefi °f complex numbers such that 
for all k £ N, A nfc G E rik and lim X rik = A. 

k— >oo 

• The set lim E n (limit inferior) is the set of all complex numbers A G C such that there exists 
a sequence (\ n ) n <£N of complex numbers such that for all n £ N, X n G E n and lim A„ = A. 

• If lim E n = lim E n , then lim E n := lim E n = lim E n . 

n-»oo n->-oo n->oo n->-oo n->oo 

2.2 Consistent and non-consistent approximations 

Definition 2.2. An approximation (7^) n eN of a self-adjoint operator A is a sequence such that, 
for all n G N, 

In ■ [X n , a^i , TT^ti ) , 

where 
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• (X n ) n en is a, sequence of finite dimensional subspaces ofQ(A); 

• («Ti)ngM is a sequence of symmetric bilinear forms on Q(A); 

• ("*n)n6N is a sequence of symmetric bilinear forms on Q(A) such that the restriction of m n 
to X n forms a scalar product on X n . We denote by \\ ■ \\x n the associated norm: V</>„ £ X n . 

\\<t>n\\Xn = m n{^n^n) 1/2 ■ 

The approximation (7^)neN is called consistent if, for any (ip, A) solution of (OJ), 

V^„ £ X n , a n (ip, (j> n ) — Xm n {i>, (f> n ), 

and non- consistent otherwise. 

The approximation (7n)«eN is referred to as a standard Galerkin method if, for all n G N, 
a n = a and m n = to. Standard Galerkin methods are obviously consistent. 

If (T n )n£N is an approximation of A, we denote by A n and M n the m-symmetric (i.e. symmetric 
w.r.t. the scalar product to(-, •)) linear operators on X n defined by: V<^ n , ip n £ X n , 

to (^4„0„, ipn) = a n (cj) n ,ip n ), 
to [M n cj) n ,ij) n ) = m n (<p n ,ip n ). 

Since to„ is a scalar product on X n , the operator A4 n is invertible and we can define the operator 

A n = M~ 1/2 A n M~ 1/2 

on X n , which is m-symmetric as well. The generalized eigenvalue problem 

J find (ip n , A n ) £ X n x R such that HV'nllx = 1 and 

\ V0„ € Xn, On(lpn, <f>n) = X n m n (lpn, 4>n), 

1/2 

is then equivalent, through the change of variable = M n ip n , to the eigenvalue problem 
f find (£„, A„) € X n x R such that ||£ n ||f/ = 1 and 

The main objective of this work is to provide sufficient conditions on such potentially non- 
consistent approximations (7^) n eN so that the discrete eigenvalues of A and the associated eigen- 
vectors are well-approximated in a certain sense by eigenvalues and eigenvectors of the discretized 
problems ([2]). We wish to provide a framework which will enable us to deal with the supercell 
method for perturbed periodic linear Schrodinger operators described in Section [SJ 



3 An abstract convergence result 
3.1 The general case 

Let us consider an approximation (Tn) n eN of A satisfying the following assumptions: 



(Al) € Q(A), 



i-n 



Q(A) 
X n 



0; 



Q(A) n-voo 

(A2) there exists < 7 < T < 00 such that for all n £ N and all ipm<t>n £ X n , 

\On(lpn,<t>n)\ < F\\lp n || q( A ) \\ <j) n \\ q( A ) ] 

(A3) for any compact subset K C C, if there exists a subsequence (T nk )k&i of [T n )nen such that 
dist (K, a(A nk )) > ax for some ax > independent of k e N, then there exists ck > such 
that for all /j, £ K and all k £ N, 



inf sup 

wn k ex nh Vnk ex„ k 



\( a n k H m n k rifc j ^nfc J 
||™tiJ|q(A)|KJ|q(A) 
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(A4) there exist k £ M+ and, for each n £N, two symmetric bilinear forms a n and m n on Q(A), 
and four seminorms r", r™, s°, s™ on such that Vc/) n ,ip n £ X n , 



y\\4>n\\u < fn n {ij) n ,ij) n ) < riiV'nllli) 

|an(-0n,0n)| < r|| ^Vi II Q(A) II <An 1 1 Q(A) ) 



and V<£, V S 



|(a-a„)(0,V)| <<(0KW, |(m-m„)(W)|<C(CW, 
r»(0) < «||0|| q( a), C(^) < 

5(0) — ► 0, C(^) — ► 0, 



r 



n— >oc 



and 



(a„ - a n )(n^ A V)W„)| ^ a ^ iN |(m„ -m„)(II^ A) ^,u; n )| 



sup - — ,| <<w, sup - — <<m 

<(</>) — ► o, — ► o. 

n— ^oo n— >-oo 



Before stating our main result, let us comment on these assumptions. 

Conditions (Al) and (A2) are classical. The former means that any ip £ Q(A) can be approxi- 
mated in Q(A) by a sequence (V'n)neN such that i/j n £ X n for each n £ N. The latter ensures that, 
uniformly in n, the norms || • \\x n and || • \\n are equivalent on X n , and the bilinear forms a n are 
continuous on X n , the space X n being endowed with the norm || • || q(a) ■ 

Assumption (A3) is important in our proof since it enables us to apply Strang's lemma (see 
Section [8]) with a uniform discrete inf-sup condition. 

For the supercell approximation, we will prove a stronger result: 

(A3') for any compact subset K C C, there exists ck > such that for all n £ N and all [i £ K , 

. \(a n - um n )(w n ,v n )\ . 

mf sup — 1, 1, 1, > c K min(l,dist(/x,CT(A„))). 

w n eX nVnEXn \\Wn\\Q(A)\\Vn\\Q(A) 



It is easily checked that (A3') implies (A3). 

Let us finally comment on condition (A4) in the perspective of the analysis of the supercell 
method with numerical integration addressed in Section [3] In the latter setting, the introduction 
of the bilinear forms a n and fh n aims at separating in the error bounds of Theorem 13.11 the 
contributions inherently due to the supercell method (truncation of the domain and artificial 
periodic boundary conditions) and those due to numerical integration. We postpone until Section[5] 
the precise definitions of a n , m„, a n and m n in this context. 

Note that (A4) implies that the approximation (7" n ) n eN is weakly consistent in the sense that 
for all <j> £ Q(A), the consistency errors r%(cf)), r™(cf)), s^(</>) and s™(0) converge to as n goes to 
infinity. 

We are now in position to state our main result. 

Theorem 3.1. Let A be a self-adjoint operator onH, X £ cr^A) a discrete eigenvalue of A with 
multiplicity q, and {7~ n )neN an- approximation of A satisfying assumptions (Al)-(A^). Then, 

1. Convergence of the eigenvalues 

A G hmcr(A„). (3) 

2. A priori error estimates in the absence of spectral pollution 
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Assume that 



(Bl) 3s > s.t. (A - e, A + e) n a (A) = {A} and lim a(A n ) n (A - e, A + e) = {A}. 

n— ►oo 

Le£ P := l^j.^) be the orthogonal projection on Ker(^4 — A) and 

V n := ix„Ai- 1/2 l (i _ E/2 a +£/2) (4)^y 2 ^, 

Rank("P„) > q, 
and there exists C € R+ smc/i iftai, /or n large enough, 



\\{V-V n )V\\c(nMA))<C 



1 - U Q X [ A) ) V 



C{H,Q(A)) 



(4) 
(5) 



with 





:= sup 

T/jeRan(T>), ||i/>||w 


=1 




:= sup 

^SRan(P), HV-II-H 


=1 




:= sup 

■4>eRnn(V), ||i/>||w 


=1 




:= sup 

T/jeRan(T>), ||-</>||w 


=1 



If we assume in addition that 

(B2) for n large enough, Rank("P n ) = q, 
then there exists C £ K+ such that, for n large enough, 



< C 



i - n^ A) ) v 



{V - Vn)Vn\\c(-H,Q(A)) 

|A„-A|<cf(||(i-n^>) 



max 

A„eo-(A„)n(A-e/2 : A+e/2) 



It is easy to check that P n := M. n l(A-e/2,A+e/2) (A n )A4n* is the m„-orthogonal projection 
of X n onto the space Y n C X„ spanned by the eigenvectors of ((2J associated with the eigenvalues 
belonging to the interval (A — e/2, A + e/2). The operator V n = ix n Pni*x € i s therefore a 

(non-orthogonal) projection on the finite dimensional space ix n Y n C "H. 

Theorcm l3 . 1 1 implies that, if (7^) n eN is an approximation of the operator A satisfying (A1)-(A4), 
for all discrete eigenvalue A of A, there exists a sequence (X n )neN of elements of a(A n ) converging 
to A. Assumption (Bl) states that there is no spurious eigenvalues in the vicinity of A. Estimate 
([5|) shows that under assumption (Bl), for each eigenvector ip of A associated with the discrete 
eigenvalue A, there exists a sequence (tp n )neN of elements of Ran("P„) which strongly converges 
towards -0 in Q(A). 

On the other hand, there may a priori exist a sequence (ipn)neN of normalized elements of 
Ran('P n ) weakly converging in H towards a vector that is not an eigenvector of A associated 
with A. This is excluded when we make the additional assumption (B2). Assumption (B2) means 
that, for n large enough, the sum of the multiplicities of the eigenvalues of A n close to A is equal to 
the multiplicity q of A. Under this assumption, if (tp n )neN is a sequence of vectors of % such that 
for each n large enough, tp n is an "H-normalizcd eigenvector of A n associated with an eigenvalue 
A n G (A — e, A + e), and if (ipn)n<£N weakly converges in % towards some ip € H, then estimate (j6|) 
implies that ip is a "H-normalized eigenvector of A associated with the eigenvalue A and that the 
convergence of (V'n)neN to ip holds strongly in Q(A). 



C(H,Q(A)) 

V 

C(H,Q{A)) 
,1/2 



n a n + K™+S a n + S™\, (6) 



(7) 



6 



Lastly, estimate ([7]) shows that when a n = a n and fh n — m n (which is the case in the supercell 
model when numerical integration errors are neglected), then 5^ = S™ = 0, and the convergence 
rate of the eigenvalues is twice the convergence rate of the eigenvectors measured in the Q(A) 
norm. Such a doubling of the convergence rate is expected in variational approximations of linear 
eigenvalue problems (see e.g. [5]). 



3.2 Standard Galerkin method 

Let us now consider the special case when for all n £ N, T n = (X n ,a,m) where (X n ) ne ^ is a 
sequence of finite dimensional subspaces of Q(A) satisfying (Al). In this case, for all n £ N, 
A n = A\x n , M-n is the identity operator, and 

Vn = ix n ^(X-e/2,\+e/2){A\x n )i*x n 

is an orthogonal projector with respect to the scalar product m. 

In this setting, (A2) and (A4) are obviously satisfied, and (A3) and (A3') respectively read 

(C3) for any compact subset K C C, if there exists a subsequence (X nk )k e m of (X n ) ne ^ such that 
dist (K, a(A\x n )) > Q-k for some ax > independent of k £ N, then there exists ck > 
such that for all n £ K and all k £ N, 

|(a-HK,^)| . 

ml sup > Ck, 

Wn h eX nh v nk eX„ k \ \ W n J | Q( A ) \ \ «n fc 1 1 Q ( A) 

and 

(C3') for all compact subset K C C, there exists Ck > such that for all n £ N and all fx £ K, 
. \(a- fj,m)(w n ,v n )\ . 

mf sup T i, 1^ — I, > c K mm(l,<bst(ii,cr(A\x n ))). 

w n ex nVneXn \\w n \\ Q{A) \\v n \\ Q{A) 

It is proved in [T5[ that, when A is semibounded, (C3'), and thus (C3), automatically hold. On the 
other hand, when A is not semibounded, (C3) is not always satisfied. An explicit counterexample 
is given in [15] . 

The formulation of Theorem 13.11 simplifies in this case as follows: 

Corollary 3.1. Let A be a self-adjoint operator onH, X £ Ud{A) a discrete eigenvalue of A with 
multiplicity q, and (A„)„ s n a, sequence of finite dimensional subspaces of Q(A) such that 



VV> e Q(A), 



— > 0. 

Q(A) ra->oo 



Let us assume that either A is semibounded or (A„)„ e N satisfies assumption (C3). Then, 

1. Convergence of the eigenvalues 

A e limo-(A|y n ). 

2. A priori error estimates in the absence of spectral pollution 

Assume that 

(Dl) 3s > s.t. (A -e, A + e) n a(A) = {A} and Urn a(A n ) n (A - e, A + e) = {A}. 

n— f oo 

LetV :— be the orthogonal projection onKei(A-X) andV n '■= ix n ^(\-e/2.\+e/2){A\x n )i*x 

Then, 

Rank(7 : ' n ) > q, 
and there exists C £ R + such that, for n large enough, 

ll(^-WII^))<c||(i-nr)^| £( ^ Q(A)) - (8) 



7 



If we assume in addition that 



(D2) for n large enough, Rankin) = q 
then there exists C S K+ such that, for n large enough, 

\\{V-V n )V n \\ C{ U,Q(A))<C 



I- Tiff) V 



max 

A„SCT(A„)n(A-e/2,A+e/2) 



\x n -x\<c(\\(i-iiff)v 



C{H,Q(A)) 



(9) 
(10) 



The estimates ©, © and (ITU1) are optimal. They are similar to the ones proved in [TH 120"! |2"T] . 
but our assumptions on the sequence of discretized operators A\x n are different. In [14], these 
estimates are proved under the condition 



S(A,A\ 



where 



S(A,A\ X J 

In [20], the assumptions are that A is invertible and 



sup inf 

£D(A), \\<f>\\n + \\A4>\\n = l' l '™ eX ™ 



\A~ 



sup inf 

i„ex„ w n ex n 



0. 



U-<pn\\n + \\A<f>-A\ Xn (f>„\\H- 



\Q(A) 



(11) 



IKIIq(A) 



(12) 



Each of the conditions (fTT|) and (fl"2|) ensures that (Dl) and (D2) hold for any discrete eigenvalue 
of A. In the case when A is semibounded, (C3) is automatically satisfied, so that our assumptions 
boil down to (Al), (Dl) and (D2). These three conditions are weaker than those in [141 EH] , 
and more easy to check in some settings, as will be seen in Section[5]on the example of the supercell 
method. On the other hand, when A is not semibounded, the precise relationship between condition 
(C3) and (fTTjl and (fT2|) is still unclear to us. 



4 Proof of Theorem 13.11 
4.1 Proof of (SD 

Let us argue by contradiction and assume that there exists a subsequence (7~ nk )k<£N and r\ > such 
that (A - 77, A + r\) n a (A) = {A} and 

VfceN, dist{\, a(A nk )) >r). (13) 

Let ip £ D(A) be a ^-normalized eigenvector of A associated with the discrete eigenvalue A and 
[i := A + tj/2. As (fi - §, [i + §) H a (A) = 0, it holds 

\v — fi\ 

a := mm — - > 0. 

uea(A) 1 + \u\ 

Let us consider the auxiliary problem 

J find u e Q(A) such that , , 

[ Vu <E Q(A), (a — /jm)(u, v) = (A — n)m(Tp, v). 

The bilinear form a — /im is continuous on Q(A) and satisfies \\a — ^Tn\\c(Q(A)y.Q(A)) < 1 + | /xj - 
The linear form / : Q{A) 3 v <— > (A — /i)m(?/>, u) is also continuous. Furthermore, as /i ^ c(^4), if 
w G Q(^4) is such that (a — fim)(v, w) — for all w € Q(^4), then necessarily v — 0. Lastly, 

. ,• l(a — nm){v, w)\ . \v — a\ 

inf sup V „ P ./y 71 > min I ^ = a. 

u,GQ(A)„ eg(A) |M|q(A)IM|q(A) vecr(A) 1 + |z/| 
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Thus, applying Banach-Necas-Babuska's theorem (see Section [8]), problem (fT4)l is well-posed. 
Clearly, its unique solution is u = tjj. 

Let us now introduce the following sequence of discretized problems for k G N: 



find u rik <E X nk such that 

Vw„ fc S X nk , (a nk - jUm„J(u nfc ,w n J = (A - fi)m nk (ll^ip,v nk ^ . 



(15) 



From (fTB")) and assumption (A3) (for if = {/i} and = 77/2), we deduce the discrete inf-sup 
condition 

\{a nk - ^m nk )(v nk ,w nk )\ 



\/k € N, inf sup 



> c> 0. 



\\ v n k \\Q(A)\\Wn k \\Q(A) 

Thus, by Strang's lemma (see Section [5]) and assumptions (A2), (A3) and (A4), for all k e N, 



\\ip - u 



n k \\Q(A) 



< 



>1 



sup 



2c t ,„ fcS x„ fc \\ v n k \\Q(A) 



inf 

w nu eX. 



n , .I 1 («» fc ~ g + ~ m nk ) w» fc , v nk ) 

\\1p-Wn h \\Q{A) + - Slip - - 

°v nk eX n . \\Vn k \\Q(A) 



< 



2c 



n 



c+i + ImI 



Q(A) 



C* 



+ C (V) + C k (V) + a* (V) + sZ ty)) , 



Q(A) 



where C € K+ is a constant independent of k. The above inequality implies that the sequence 
{u nk )keN strongly converges to ip in Q(A), from which we infer that 



, lim Pxl^V - U nk \\ Q (A) = 0. 

On the other hand, (|15p yields 

Vu„ fe £ A Ilfc , (a„ fc - Am„J(u„ fe ,u n J = (A - |ii)m nj , (n^V 
The above equality also reads 

(A nk - \){M]/ k 2 u nk ) - (A - ^A^f (n^V - «n, 
It then follows from (Al), (A2) and ([16]) that 

lim || - A)(My>„J||« = and lim inf ||A< 2 U „J|^ > 7 > 0, 

fc— >oo fc— s-oo 



(16) 



which proves that dist(A, cr{A nk )) — >■ and contradicts (fT3|) . 
4.2 Proof of dU) and (JSJ 

By assumption (Bl), the approximation (7^) n eN is such that lim <r(A n ) H (A — e, A + e) = {A}. 
Hence, for n large enough, 

cr(^ n ) n ((A - 2e/3, A - e/3) U (A + e/3, A + 2e/3)) = 0, 

so that the circle C in the complex plane centered at A and of radius e/2 is such that dist(C, cr(A„)) > 
e/6. This implies in particular that, for n large enough, 
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Consequently, for all tp £ Ran("P), it holds 

(V - V n )i> = ^ &> ([z - Ay V - ix n M- x l\z - A n y l M]! 2 i* x Jj) 



dz. 



In the following, C will denote a constant independent of n £ N* and z £ C, which may change 
along the calculations. 

For z £ C, we introduce the auxiliary problem 

find u z £ Q{A) such that 
Vu £ Q{A), (zm — a){u z , v) ~ m (ip, v) , 

whose unique solution is u z — (z — A)^ 1 ^ ~ since tp £ Ran('P). We also introduce the 

discretized problem 

find £ X n such that 

Vw„ £ X n , (zm„ - a„)«, v n ) = m n (n^V; «n) - 

whose unique solution is u z n — ixnMn. 1 ^' '{z — An)^ 1 Ail/ 2 i* x ^tp . From assumption (A3), since C is 
a compact subset of C, there exists c > such that for all z £ C and n £ N, 

. |(a„ - zm„)(v„,u; n )| 

mf sup — — - > c. 

w n ex n VneXn \\v n \\ QiA) \\w n \\ Q{A) 

Reasoning as in Section B~Tl we infer from Strang's lemma, assumptions (A2)-(A4) and the fact 
that r°, r™, s° and are semi-norms, that for all z £ C, 

„ z zu . 1 |m(V',w„) - m„(n^ V»«n)| 
W — "n < - sup j| - 

c »„EX„ IK Q(A) 



inf 

io„GJf„ 



c + l + \z\.. , .. 1 

II" - w n\\Q(A) + - SU P 

C v n €X„ 



hi - a) + z(m n - m)](w n ,v n )\ 



c C »„€X„ |KIIq(A) 

< ^C(V) + ~ || V - n« v|| w + + ~ |n^ A) v - 

c+ 1 + |z| 

U - lly U 

0(A) 



< 



+i (< (n<£V) + + \z\nr™ (n£V) + |z|C(u z ) 
^) + - *-i& A) 1> + c + (1 + N)(1 + k2) 

C " Vr/ C X " ^ Q(A) C 

+- + + MO 2 ) + |*|<(«*)) 



0(A) 



c* 



v ' 0(A) 



since u 



z-X 



. Thus, for all z £ C, 



(z - - ix n M~ 1/2 {z - A„)- 1 7Wy 2 ^„V 



Q(A) 



< c 



l-n 



0(A) 



Q(A) 



Since C is of finite length, we obtain that, for n large enough, for all tp £ Ran('P), 



\\(P- V n )tP\\ Q{A) <C 



i-n^ A >) 



tP 



Q(A) 



which readily leads to 

Let us finally consider a "H-orthonormal basis • • • , £ g ) of Ran(7') = Ker(A — A). Since for 
all 1 < i < q, V n d — > VCi = d strongly in H, the family (V n Ci, • • • )V n Cq) is free for n large 

n— »oo 

enough, so that Rank("P„) > q. 
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4.3 Proof of (ED and flJD 

We just have shown that for all 1 < i < q, V n C,i — > = (i strongly in H. Under the additional 

n— f oo 

assumption that, for n large enough, Rank(7 : '„) = q, this implies that there exists no G N, such 
that, for n > rig, (PnCij ' ' ' ^nCq) forms a basis of Ran(V n ), with 

3 1 

mm ||PnCi|& > T and max \m (V n &, 0)1 < 7- 

Thus, any £ n € Ran("P„) can be decomposed as 

i 

£n = $>*(&.)P»Ci, 

i=l 

the coefficients (a 1 (^„), • • • , of in the basis (VnCi, ■ ■ ■ jVnCq) being such that 

max |a l (C„)| < 2||&,||«. 

l<i<q 

We have 

q I q 



i=l 
q 



i=l 



and we deduce from ([5]) that for all 1 < i < q, 

IICi-^Cillg W <c(||(i-n^)p 

Hence, 

V£„ € Ran(7>„), ||7>£ r -CnllQ(A) < C 



C(n,Q(A)) 



C(n,Q(A)) 



where the constant C is independent of n. Besides, it also follows from (A2) and the definition of 
V n that 



VneN, ||P„|U W < W- 



Therefore, 

\\{V-V n )V n \\ m ,Q{A)) < 

< c 

and © is proved. 



Il^n - £«||q(A) ,, 

SUP FnlU(«) 



(i-n£ A) ) 



For each n large enough, let (^„, A„) € X n x R be a solution to the generalized eigenvalue 
problem © such that A„ € (A - e/2, A + e/2), <A„ = and Xn = = yffifa- It 



follows from ([6]) that 

ll^n-^nlk < WPK-<Pn\\ Q( A) < C (J\ (l - n^ A) ) 



£(«,Q(A)) 



0, 
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from which we infer that H'P^nllw — > 1, (<An.)«.eN is bounded in Q(A), \\<fi n — Xn||Q(.A) ~> 0, and 



llXn - <Mq(A) 



< 



||^. 



n\\H 



-Vcj> r 



Q{A) 



+ \\P4>n ~ 4>n\\Q(A) 



< ll^ ,< / > n|l« 1 |l < /'™ - 'P<t> n \\H\\'P<t>n\\Q{A) + WP<Pn ~ <Mq(A) 



c 



(i-n 



Q(A)\ 
x„ ) 



V 



C(U,Q(A)) 



Besides, it holds 



|A„-A| = |a n (V>n,^n) - a(Xn,Xn)\ 

a n ((f) n ,(f> n ) 

< : — — a{<Pn,(Pn) 



+ \a((p n ,(f> n ) - a(Xn,Xn)\- 



m n ((j)n, 4>n) 

On the one hand, we have 

|<l(0nj0n) - a(Xn,Xn)\ = \a{4>n ~ Xn, <Pn ~ Xn) + 2a(x„, <t>n - Xn) I 

= |a(0n - Xn, 4>n - Xn) + 2Xm(Xn, 4>n ~ Xn)\ 
= \a{4>n - XnAn ~ Xn) ~ A||Xn - 0n\\n\ 
< C||0n-Xn||Q(A)- 

On the other hand, 

\(a — a„)(0„, (j> n )\ < | (a — a n ){4>n, 4>n)\ + |( )(<^i,<MI 

< ^n(<^n) 2 + | (fin ~ a n ) - H%[ A) Xn , 0„ - Tlf A) X n) 

+2 |(a„ - a„) (n^ A) X n, - ng^x 
+ |(a„-a n )(n^x n ,n^xn) 

< «(Xn) + «||&, - XrjQ(A)) 2 + (T + AC 2 ) 



n Q{A) Y 



Q(A) 



Q(A) 



U Q{A) Y 



Q(A) 



s n(Xn) 



c 



c 



(r a n (Xn) + Un- Xn\\ Q( A) + || (l " Tiff) Xn\\] + <(X») 



i-n« A) ) 



v 



C{H,Q{A)) 



•5° 



and a similar calculation leads to 

|W„( <t>n,(t>n)- 1| = |TOri(^n, <An) ~ ™(0n, 0n)|> 



c 



i-nf A) ) 



C{H,Q{A)) 



Consequently, 

a n (<j>n,(t>n) 



m n (4> n ,4> n ) 



a{4>n,<t>n) 



, I (a - a n )(</>«, </>n) I , , 
- 77 H°(0ni0rj| 



m n (4> n , 4> n 



m n (4> n , (f> n ) 



< 7 (| (a - a n ){4>n,4>n)\ + \a{4> 

)| |m„(0 n , 4> n ) — 1|) 



< C 

Collecting the above results, we obtain 



i-nf A) ) 



v 



C(H,Q(A)) 



|A-A„| <C 
which proves estimate ([T]>. 



i-n« (A) ) 



A'„ 



C(U,Q(A)) 



111 +11 



n 1 n 



■5° +S, 



n 1 n 
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5 Application to the supercell method 



The aim of this section is to show that the theoretical framework presented in Section [3] can 
be applied to the numerical analysis of the supercell method for perturbed periodic Schrodinger 
operators. 

Note that the supercell method was previously studied from a mathematical viewpoint by 
Soussi [26] , for the special case of a two-dimensional periodic Schrodinger operator in the presence 
of a compactly supported perturbation W of the form W{x) = w\q,{x), where w is a real constant 
and Q, a bounded domain of M 2 . 

5.1 The supercell method with exact integration 

Let TZ be a periodic lattice of M. d , TZ* its reciprocal lattice and T a unit cell of TZ such that is in the 
interior of T. Typically, in the case of the cubic lattice TZ = Z d , TZ* = 2n1 d and T = (-1/2, l/2] d 
is an admissible unit cell. 

Let us introduce the perturbed periodic Schrodinger operator 

A := -A + Vp er + W, 

where A is the Laplace operator, V pel a real- valued 7^-periodic function of .Lf oc (R d ), with p = 2 if 
d<3,p>2iftf = 4 and j3 = ff/2ford>5, and W £ L°°(R d ) a real- valued function such that 
W(x) —> 0. 

foo 

The operator A is self-adjoint and bounded from below on % ■— L 2 (R d ), endowed with its 
natural inner product 

V</>, € Ti, m(4>,ip) := / <fnp, 
JR d 

with domain D(A) = H 2 (M, d ) and form domain Q(A) = i/ 1 (K d ). The associated bilinear form 
a(-, •) is defined by 

V0,^eQ(A), a(<f>,il>):= f V<£-VV> + / (F per + W)H). 

JR d JR d 

We denote by A := — A + V per the corresponding periodic Schrodinger operator on L 2 (M. d ). 

The supercell method is the current state-of-the-art technique in solid state physics to compute 
the spectrum of the operator A. For L £ N*, we denote by Tl := LT the supercell of size L and 
set 

L 2 pcl {T L ) := {u L G Ll c (R d ) \ u L Lft-pcriodic} , 

ffperOi) := {u L e i^ r (ri) | V Ui G (L* or (r L )) d } , 

Cper(ri) := {« L G C°(M d ) | u L Lft-periodic} , 
LpevpL) ■= {u L e L°°(R d )\u L LTZ-periodic} . 
For u L G L 2 C1 .(T L ) and fe G L^ 1 TZ* , we denote by 

«i(*0 : = ^^1/2 ^ u L {x)e lk - x dx 

the Fourier coefficient of corresponding to the fc mode. For r G R, the Sobolev space H^^l) 
can be defined as 

H r pci (r L ) :=\u L e L 2 pa .(r L ) | £ (i + |fc| 2 ) r \u L (k)\ 2 <00 \. 

The supercell method relies on the resolution of the following (non-consistent and non-conforming) 
eigenvalue problem: 

J find (ul,n,^l,n) € ii.iV x K such that 

\ Vt>£,jv G F^jv, a L (u LyN ,v LtN ) = X L ,Nm L {u l ,n,Vl,n) > 
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where 

Vu L ,v L £ Lp Cr (r L ), fh L (u L ,v L ) := / u^Ul, 



Vu L ,v L £ iJp Cr (r L ), a L (u L ,v L ) ■■= / Vml ■ Vul + / (Vp Cr + W)ulvl, 

and Yj,^ is a finite dimensional subspace of iJp C1 .(r^). 

We set i/^jv = Hl\y L}N , where Hl denotes the unique self-adjoint operator on Lp er (rL) 
associated with the quadratic form cll- We have D{Hl) = H 2 er (TL) and 

Vu L £ -ffp Cr (r L ), HlUl = ~Au L + (Vper + W L )w L , 

where W^l € Lp^fTi) denotes the L7\L-periodic extension of W|n- 

For the sake of clarity, our analysis will be restricted to the case of the cubic lattice 1Z = Z d 
and the planewave discretization method, for which 

Yl,n ■= I 51 Cke-L,k | Vfc, c_ fc = < 

[fc6L-i7?.* | \k\K2nNL- 1 

where eL,k{x) := |ri| _1 / 2 e lfc :c . We denote by IIe,^ the orthogonal projection of L^^Tl) on Yn,l 
for the Lp er (ri) inner product (actually Ul.n is also the orthogonal projection of ffp er (rj,) on 
Yjv,i for the i?p Cr (r i ) inner product, for any set). 

The discretization spaces Yr, jv possess the following properties: 

Vul,7V G Xl,AT, IIl,7V (-Aui,iv) = -AuL t ff, 

and for all real numbers r and s such that < r < s, there exists a constant C £ R+ such that for 
all LeW and all u L £ H* m (r L ), 



\\U L - H L ,NUL\\H; m (rt) <C[jf) \\ u L\\H^(r L )- ( 17 ) 

As in [8], we will assume that V pol belongs to the functional space Z per (T) (denoted by A^ per (r) 
in [8]), defined by 

2 por (r) := \V £ LpUD I \\V\\z p „ { r) ■= sup sup ^ WL ^U^) < +0O | 

The space 2 per (r) is a normed space and the space of the 7^-periodic functions of class C°° is 
dense in Z por (r). 

Our main result concerning the supercell method in the absence of numerical integration error 
is the following: 

Theorem 5.1. Assume that V peI € Z pQr (T) and that W £ L°°(R d ) with W(x) — > 0. Let 

\x\— >oo 

N L 

(-^l)lsn* be a sequence of integers such that —— — > +oo. Then, 

L L— >oo 



1. Absence of pollution 



2. A priori error estimates 



Km a(H L , Nti ) = a(A). (18) 

L— f oo 



Assume that, in addition, V pcl £ Hp CI 2 (T) and W £ H r 2 (R d ), for some r > 2. Let A be a discrete 
eigenvalue of A ande > be such that a(A)f] (A-e, A + e) = {A}. Let V := 1{ X }(A) be the L 2 (R d )- 
orthogonal projection onto the eigenspace of A associated with A and^L := l(A-e/2,A+e/2) {Hl,n l ) 
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the Lp er (r l)- orthogonal spectral projection of Hi / ^n l associated with the eigenvalues belonging to 
the interval (A — e/2, A + s/2). Consider finally a sequence of cut-off functions (xl)lgn* such that 

< XL < 1 on R d , xl = 1 on T L , Su PP ( X l) C (L + VI)r, ||Vxl||l~ < c, (19) 

for some constant c G R+ independent o/ieN*. 

Then, Ran('P) C if r (K d ), and there exists C, 8 > suc/i i/iai /or L large enough, 



Tr(V) = Tr(q3 L ), (20) 
sup inl^ \\^- XLUL \\ Him <C\e- SL +(^-) j, (21) 



^GRanf-p), || V|| L 2 {Rd) =1 Mi eRan^i ) 



L ^ r - lN 



sup , inf ||V - Xi^||H 1M < C e" 6L + ( T7- ) I , (22) 

max |Al — A| < C ( e~ 5L + ( -rf-\ ] ■ (23) 

5.2 The supercell method with numerical integration 

In general, the computation of the integral J rL {V pel + W)ulvl with ul,vl G Yl_n l cannot be 
carried out explicitly, and a numerical integration procedure is needed. We assume in this section 
that V pcl and W are continuous functions. 

For Mef and ml G Cp CT (rL), we denote by u L ' the discrete Fourier transform of ul on 
the cartesian grid Ql,m '■= -jj^ a 1 ■ Recall that if 

ul — ^ u L (k)eL,k, 

the discrete Fourier transform of mj, is the ML _1 7\L*-pcriodic sequence u FFT,M = {u FFT,M ' (k)j 
where 

u F l FT ' M W = W E u L {x)e~^ = \T L \~^ Y, Mk + MK). 

xeg L ,Mnr L KeL-^n* 

We now introduce the subspaces 



wl D M 



Span \e lty \ I G 27rL _1 Z ) |Z| < ^ ^ } (M odd) 

Span je^ | Z G 2ttL- 1 Z, |Z| < ^ ( ) ) ® C (e™ Mv l L + e""^^) (M even 



and denote by Wl,m the d-tensor product space Wl,m ■— Wl° M ®- ■ - ®Wl D M . In particular, when 
M is odd, 

Wlm = Span! e L , fc , k G i" 1 ^*, |fc|oo < 271-iT 1 



It is then possible to define the interpolation projector Tl,m from Cp 0r (r i ) onto VFl.m by [Zl,m( u l)} 
Ul(x) for all x G Gl,m- ^ n particular, when M is odd, we have the simple relation 

Xl, M ^l) = \Tl\ 1/2 E u F L FT - M (k)e L , k . 

keL-in* | |fc| DO <2 7r L- 1 (i^) 

It is easy to check that if the function mj, is real-valued, then so is the function 1l,m(ul)- 
Besides, when M > AN + 1, it holds that for all ul 7 vl G Yl,n, 



/ Tl,m(V l u l v l ) = / Il.m (Vl)u l v l , 

JTl JTl 
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for any V L G L 2 el (T L ). 

The supercell method with numerical integration then consists in considering the following 
eigenvalue problem for a given M > 4N + 1 , 

find (ul,n, Al,jv) € Yl.n x K such that 
Vvl,n € yt.jv, o,l,m(ul,n,vl,n) = >^L,Nm L (ul,n, vl,n) , 

where 

yu L ,v L G ffp Cr (r_L), a L ,M{uL,v L ) ■= I Vu L -Vv L + I lL,M(V peT + Wl)u l v l , 

and where Wl is the L7\L-periodic extension of £l W\r L , £l being a C[' -i](K d ) cut-off function such 
that < £l < 1, £l = 1 on r L -i, Supp(^) C (L - 1/2) T, and the sequences (\\d a Len , 
are uniformly bounded in L, for all \a\ < [r — 1] (here and above, [r — 1] denotes the integer part 
ofr-1). 

As in the preceding section, we denote by Hl,n,m = Hl,m\y l n , where Hl,m is the unique 
self-adjoint operator on Lp er (ri) with domain D(Hl,m) — H 2 ci: (Tl) associated with the quadratic 
form cll,m- 

N L 

Theorem 5.2. Let (Nl)l£j$* and (Gl)lgn* be sequences of integers such that — > +oo and 

G L — > +oo, and M L := LG L . We assume that V pcl G C° (r) n H;~ 2 (T) and W G C°(R d ) n 
ff r - 2 (R d ) for some r > 2. TTien, 

i. Absence of pollution 



lim a(H LtNLt M L ) =tr(A). (24) 



2. ^4 priori error estimates 



Let X be a discrete eigenvalue of A and e > be such that a(A) n (A — e, A + e) = {A}. Le£ 
■p := lr^i(A) be the L 2 (M. d )- orthogonal spectral projection onto the eigenspace of A associated with 
A, and^ L := l(\- £ /2,x+e/2){H L ,N L ,M L ) the L 2 cr (r L )- orthogonal spectral projection of H l . Nl . Ml 
associated with the eigenvalues belonging to the interval (A — e/2, A + e/2). We finally consider a 
sequence (xl)lgn* of cut-off functions such that 

< XL < 1 on R d , xl = 1 on T L , Su PP ( X l) C (L + VZ)T, ||Vxi|U~ < c, 

for some constant c G R + independent of L G N* . 

TTien, Ran('P) C iJ r (R d ), and i/iere exists C,8 > such that for L large enough, 

Tr(V) = Tr(<p L ), (25) 

sup inf \\i>-XLUL\\w<m*)<C(ei{L) + e 2 (L)), (26) 

V-GRanCP), ||^|| i 2 (s d ) =l u ^ eRan (^) 

sup inf ||^-xi,UL||jfi(Rd) <C(ei(L)+e 2 (L)), (27) 

UiSRanfqJi), UutU^a (r , , =1 */>eRan(P) 

max |Ai-A|<C(ei(L) 2 + e 2 (L)), (28) 

w/iere 

eiW == e"« + ( j^)' 1 e a (L) := (A)' ' + | W | (e"« + (A)") . 
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5.3 Formulation in terms of non-consistent approximations 

The supercell method can be rewritten as a non-consistent approximation of the operator A (in the 
sense introduced in Section [2.2l) . based on the approximation spaces (Xl)l£n* and the symmetric 
bilinear forms (cil)lgn* 7 (ol)l<en*, and (tol)l<en* defined for all L G N* by 



X L ■■= {xlu Li u L G Y LiNh } C ff 1 (M d ), 



and 

V0,V G ffW) 





:= 


[ V0- 

/rx, 




Sl(< 


:= 






mi( 




[ #, 





where we recall that (xl)-Lgn* is a sequence of cut-off functions satisfying (|19p. It is easily checked 
that for all L G N*, mj / (-, •) defines a scalar product on Xl- 

Let us introduce, for each L G N* , the unitary operator 

h ■ ( Y L,N L ,{-r}Ll et (rL)) (X L ,m L (-,-)) , 

ul Xlu L - 

Its adjoint (and inverse) j* L is given by: V</>l G Xl, Jl^l) = ul where ul is the L7?.-periodic 
extension of 4>l\v l - The supercell problems 

find (\l,ul) G K x Y l ,n l such that \\u L \\Ll al (r L ) = h 
V«l G yL,JV L , a L (u L ,v L ) = X L m L {uL,VL)i 

and 

find (Al,u l ) elx F Li jv l such that HwilU^^rt) = 1, 
Vv L G aL,M L (u Ll v L ) = X L m L {uL,v L ), 

are then respectively equivalent, through the change of variable ipL = 3l u Li to the generalized 
eigenproblems 

find (Al,-0l) G M x such that tol^LjV'l) — 1 and 
V0l G X L , o,l(iPl, 4>l) = ^Lm L (ipL,4>L), 

and 

find (Xl,iPl) G R x Xl such that TOl(V>z,,V>l) = 1 and 
V0l G Xl, ol^l, 0l) = AimL^L,^). 

Thus, considering the supercell method with exact and numerical integrations is equivalent to 
considering the non-consistent but conforming approximations (7l)lgn* and (7l)lgn* respectively 
defined by 

Tl = (X Ll aL,m L ) and 7l = (Xl, a L , m L )- 

Taking the same notation as in Section[3J it holds that Al = 3lHl.n l 3l an< ^ ^-l = 3lHl,n l ,m l 3l 
so that <t(A l ) = <t(Hl,Nl), = v{H l ,n l ,m l ) and, in both cases, V L = ix L iL^Li* L i*x L - Tne 

following section is devoted to the proof of Theorems 15.11 and 15.21 which are in fact corollaries of 
Theorem 13.11 We^will first check that all the assumptions of Theorem 13.11 are satisfied for the 
approximations (7l)l<en* and (7l)l£N*, and then derive more explicit expressions of the right 
hand sides of (O, ^ and ([7]) in terms of L, Nl and Ml- 

We prove in Section [6. II that the supercell method with exact integration satisfies assumptions 
(A1)-(A4). In Section lrJT2l we prove (TT5|) and (|2U1) . which imply that this method also satisfies 
assumptions (Bl) and (B2) for any discrete eigenvalue A of the operator A. Estimating the terms 
involved in estimates (JSJ, (O and ([7]) will then lead to estimates (|21[) . (|22l) and (l2lil) and conclude 
the proof of Theorem 15.11 Section 16.41 is devoted to the proof of Theorem 15.21 in which numerical 
integration errors are taken into account. 
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6 Proof of Theorem 15.11 and Theorem 15.21 



In the sequel, C will denote an arbitrary constant independent on L G N* which may vary along 
the calculations. 

6.1 Proof of (A1)-(A4) for f L = (X L ,a L ,m L ) 
Proof of (Al): Let us prove that 

V0 G H\R d ), inf ||0 - <j> L \\ Him — ► 0. 

Let <j) G H l {R d ) and e > 0. Since C c °°(]R d ) is dense in H 1 ^*), there exists 77 G C£°(K d ) such that 
H0 - ?7||fl-i(Rd) < e. Let L G N* be such that Supp(ry) C (L - y/Lo)T. For all L > L , if r\ L 
denotes the L7?.-periodic extension of r)\r L , we infer from (|T71) that 

\\VL - H-L,N L VL\\m pcr (r L ) < c j^\\Vl\\h^„(Tl) = c j^\\v\\h^) l I^ o ' 

with C € R+ independent of L. Let us then consider the sequence (0l)lgn* defined as 0l := 
Xl^l,n l Vl € for all LeN*, for which 

110- ^lIIh 1 ^) - H- vWh^ri) + Wv-XlUl^lVlWh^vl")! 

< e+ \\r) L - n.L,N L VL\\Hi„(TL) + \\XL^L,N L VL\\ H i {iL+ ^z } r\r L )- 



Furthermore, since < \l < 1, and r\jj — on (L + VZ)T \ Fl, it holds 

\\XL^LM L VL\\ 2 H1{{L+VT)rxrL) < WxiMl,n l Vl\\ 2 L 2 ((l+ ^z )V \ Tl) +^\\^Xl^l,n l Vl\\ 2 L 2 ({l+ ^z )t \ Tl) 

+ 2 \\XL^(nL,N L VL)\\ 2 L 2 ((L+ ^Z )r \ rL) 

< \\TIl,n l Vl - m\\ 2 L 2 [(L+ ^z )TXVL) + \\^{TIl,n l Vl) - ^VL\\ 2 L 2 {{L+ ^z }r \ rL) 
+ l|VxL|| L oo (R d ) ||n LiWL 77 L - 7?L|l22(( L + v ^L)r\r I/ ) 

< 3 d (4+||V X LL-( Rti ))||n iiWi 7 7i -r; L ||^ cr(rt) ^-+0. 

Hence the result. 

Proof of (A2): Let 4>l,^l G X l , and u L ,v L € Y l , Nl C Hp C1 .(T L ) such that L = xlul and 
= Xl«l- It holds 



'r 3L 



Therefore, 

Besides, 

|aL(0L,-0i)| 



,| 2 -/ x!K| 2 <3 d / \u L \^3 d f \M 2 <3 d [ I0l| 2 - 
»'r3 I/ Vrx ^rz, Jr^ 

3dll0£.|lia(R<*) < m L ((f> L ,(f> L ) < \\(/>l\\ (29) 



V0 L • V^L + / ( V pm - + W)(t> L l/jl 



< 
< 



(1 + ||VF|| L oo (R£i) )||0 L || H l (R£i) ||V'L|| ff i(R<J) + \\V VC rU L \\ L 2^ (TL) \\v L \\ L 2 cATL) 
(1 + l|W|U°°(R' i ) + \\V VC r\\z vrjI {T))\\4>L\\m(M. d )\\i } L\\m{R d )- 

Thus, assumption (A2) is satisfied. 

Proof of (A3): For all a > arbitrarily small, there exists a constant C a such that for all 
-A G iJ X (K d ), 

/ lVrll0| 2 <"/ |V0| 2 + C a / |0| 2 . (30) 

JR d J]R d JRrf 
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Besides, for all <Pl G Xl,i if <Pl = XlUl with ul G Yl j n l , it holds that 
f |V0l| 2 < 2 / _ |V X lul| 2 + |xlV Wl | 2 

iR d J(L + VZ)V 

< 2x3 d (\\\7 XL \\ Lao(md) jf M 2 +^ |Vu l | 2 

which, together with (|21)1) . yields that, for i large enough 

/ |0l| 2 +/" |V0 L | 2 <3 d+1 f / |0l| 2 +/" |V0 L | 2 V (31) 
JE d JR d \Jr L Jr L / 

Using (j3"0)l and (f3~Tj) . we obtain that for all a > arbitrarily small, there exists D a G M + such that 
for all LeW and all 4> L G 

(Upcr + VK)|(/) L | 2 < / (|U pcr | + |tU|)|0 L | 2 



< a [ \VM 2 + D a f \M 2 

JTr JT, 



This last inequality implies that there exists > independent on L G N* such that for all 

(f>L G X L , 

||^||^ (Rd) < 'i d+1 \\4>L\\ 2 HHTL) < /8(|ol(^,0 £ )| +m i (fe ) ^)). 
Thus, for all /i G C, it holds that 

|(oi - fim L )(<j) L ,ip L )\ 

mf su p in n~~n 

^eA' L ^ ie Xz, ll0i||ff 1 (R d )llw<II.F/ 1 (R d ) 

.1 . n - ^m L ){(j) L ,^ L )\ 

> — mi sup 



£ V-re^^eXx (|a L (0 L , 0l)| + m L (0 L , ^) L )) 1/2 (|S L (^ L , Vl)| + m L (^L, Vl)) 1/2 ' 
Let (Ci)i<j<(lini(x i ) ^ e an m L-orthonormal basis of Xl, such that for all 1 < i < dim(Xi), 

Hl,n l 3l<£ = 4i£Ci, 1 < * < dim(X L ), 

where {^, 1 < i < dim(Xk)} = (t(Hl_n l ). Then, any 0£ G .X^ can be expanded in the basis 
(^L)i<i<dim(X£,) : 

dim(Xt) 

4>l = c * e R ' 1 - 1 - dim (x L ), 

1=1 

and it holds that \a L (dp L , + m L (<t>L, H) < E?™^ | 2 (1 + K|). Considering 

dim(Xi) 

i=l 

we obtain that 

. |(a L - fJ.m L )((/>L, .1 . , K~mI , qo % 

ml sup j > — mf j -. (32) 

i>L£X L <p L eX L ||0L||_f/i(R^)||V'L||ffi(K ti ) P "Lea(H LiNL ) 1 + \v L \ 

Since (f3"2"|) holds for any /i G C, this implies that for any compact subset K C C, there exists a 
constant ck > such that for all LeN* and all jj, G if , 

. , |(a £ - nm L ){(j> L ,il) L )\ . 

ml sup |7— -j, :— — : > cxmm(l,dist(^,CT(iiL,JV L ))) • 

Vae^z, L eXi \\9L\\H 1 (M. d )\\' l pL\\m(K d ) 
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Thus, condition (A3'), and condition (A3), hold for the approximation (7l)lsn* 
Proof of (A4): For all <p £ ff 1 ^), we denote by 



and 



rf(d>) :-- 



rl(4>) := 



1/2 



< 



\4>\\L 2 (R d ), 



1/2 



H 2 + |V</>| 5 



< 



Then, r™ and r£ are seminorms on iJ 1 (R d ) such that for all <j) e H 1 (W i ), r™(<fr) — > and 

L—t-oc 

r l(<l>) — > 0. For all faip G tf 1 ^), it holds 



|(m - m L )(<j),il>)\ 



Let (wl)lgn* be a sequence of C°° cut-off functions such that for all L G N*, < wj, < 1, ujl = 1 
on R d \r£, Wi = on rj,_ x and the sequence (|| 'V(^L\\L ao (M. d ))Le^* 1S uniformly bounded inieN*. 
Then, for all 0, tp G i? 1 ^), 



(0-0£)(^,^)| = / V^-V^+/ (Vper + W)^ 

a d \ri, iR d \r L 



< (1+||W|U»K(0K(^)+ |y per ||a; L 



1/2 



1/2 



Using O, 



|fper||orf < § / |VK0)| 2 + C / |w L 0| a < Crl^f 



Thus, there exists k € R+ independent on L G N* such that 

|(o-o £ )(^^)|<«r£(0)r£(^). 

6.2 Absence of pollution 
Proposition 6.1. 7i /loZcfc 

cr(A) = lim <t(Hl,n l )- 

L— ¥OD 



(33) 



Besides, for any discrete eigenvalue A o/ t/ie operator A and for all e > suc/i i/iai (A — e, A + e) n 
a {A) — {A}, we have, for L large enough, 



Tr(qj L ) - Tr(7>), 
where V := 1{\}(A) and ^ L := t {X -e/2,x+e/2)( H L,N L )- 



(34) 



Let us notice that ([33]) implies that (Bl) is satisfied for any discrete eigenvalue of A, and that 
([34")) is nothing but a reformulation of (B2). We refer to [8J Theorem 3.1] for a proof of (|33j) . 

Proof of {34\l - If follows from (J33J) that (Bl) is satisfied and therefore that for n large enough, 
Tr(*PL„) > Tr(P). Let us assume that there exists an increasing sequence (ife)feeN* of integers 
such that 

Tr(y L J > q := Tr(V). 
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For all k G N, let (C^) i<i< 9 +i be an Lp Cr (rL fc )-orthonormal family of vectors of Yl^.jvz, such that 
for all 1 < i < q + 1, 

#L fc ,A^ CS = A gcg with Ag G (A — e/2, A + e/2). 



Then, for all k G N, (Xi fc Ci )i<i<g+i forms a free family of Xl h and there exists g k € Span(C)^)i<i<g+i 
such that ||<7fc||.i,2 (r Lfc ) — 1 and 

9k ■■= XL k 9k G KerfP). 
Reasoning as above, it can be easily checked that f HsfcHfl^ (r x 
that (||?fc|Ui(K d ))fc 6N 



feeN* 



is bounded, which implies 



is bounded as well. Thus, up to the extraction of a subsequence, there exists 
g G H 1 ^) n KerfP) such that g k g in H 1 {M. d ) and # fe — ► g in Lf oc (M d ). Since g k = XL k gk, 

k— foo fe— foo 

this also implies that 

g k — -> 3 strongly in if oc (M d ), 

A;— >oo 

which readily leads to 

fflz»,w £ , ~ A) g k — > -A 9 + (F per + W - A) 9 in 2?'(M d ). 
Besides, since 6 Ran(*p£ fc ) and lim cr(HL k ,N L ) = &(H), we have, 



it— »-oc 



(H Lh ,N Lk - gfc 



0. 



which, in turn, implies that 



f-ffi*,^ - A) fffc — > inZ>'(R d ). 

\ " I k—too 



Therefore, 

-Ag + (V pcl + W-\)g = 0. 

Consequently, g G Ker(P) n Ran("P) = {0}. Using similar arguments as those used in the proof 
of [5] Theorem 3.1], we infer from the fact that (g k ) k ^m strongly converges to in Lf oc (M. d ) that 

( ii~ n 9k ) is a Weyl sequence for A° = —A + Vn e r associated with A, which contradicts the 



2 (K d ) / fee 
fact that \(£a(A°) 



□ 



6.3 Proof of Theorem [5TT1 

We have proved that the supercell method with planewave discretization and exact integration 
satisfies assumptions (A1)-(A4), and that for each discrete eigenvalue located in a spectral gap of 
A, assumptions (Bl) and (B2) are satisfied. Thus, Theorem 13.11 can be applied and there exists 
C G R.+ such that for L large enough, 



\\(V-Vl)Vl\\c(l'<**),w 



TrCP) = Tr(q^), 



|| CP - 'PL)'P\\c(L 2 (R d ) ,// 1 (K d )) - C 



c 



max |Al — A| < C 

Aieo-(ffz,,jVr,)n(A-e/2,A+e/2) 



i-nf 



l-n? 



C(L 2 (R d ),H 1 (R d )) 
C(L 2 (R d ) : H 1 (^ d )) 
C(L 2 (R d ),H 1 (^ d )) 



TZ a L +11 



m 

L ) > 

2 

m 



where Vl ■= ix L iL*$Lj* L i*x L and 

K a r, 



sup r£{il)), 

sup r Z,(V0- 

^SRanCP), ||V'll_ L 2 (H d ) =l 
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Since we have 



sup inf M - XLU L \\m(R") < \\(P - T ? L)V\\c(mR<'),H 1 ( 1 st d ))> 

i/>GRan(P), 1 1 ^ 1 1 ^ 2 (IR d j = 1 "iERanOPi) 



and 



sup inf - XLU L \\m(R") < \\(P - Vl)Vl\\c(l^{ 

KiGRanOPxJ, \\u l \\ l2 =1 ^eRan(-p) 

per ^ L, > 

it just remains to prove that there exists S > independent on L such that 



i -n 



C(L 2 (R d ),H 1 (R d )) 



■R-; +-rt - c \ < •'• L ( ^ 



This estimate is based on exponential decay results for the bound states of Schrodinger opera- 
tors [25 . A real- valued function V on R d is said to lie in the class Kd if and only if 

\v(y)\ , n 
— d=2 d y = °> 



if d > 3, lim sup 

x<ER d J\x-y\<a \ x — V 

if d = 2, lim sup 



l^llndx-yr 1 ) dy = 0: 



if d = 1, 



sup 



\V(y)\dy < oo. 



zGR d J |a-y|<l 

Under our assumptions on Vp eT and W, V = V pcr + W G if^. It then follows from Theorem 
C.3.4 and Corollary C.2.3 in [25] that there exists C, 5 > such that for all L 2 (R d )-normalized 
V> G RanCP), 

VzGlT 1 , \ij{x)\<Ce~ 3Slxl and e 35M W G (L 2 (R d )) d . (35) 

For all L > 6, let ?7 L G C c °°(R d ) such that < ry L < 1, t] L = 1 on r L/2 _ 2 , Supp(?/ L ) C T L/2 -i 
and all its derivative up to the [r + l] st order are bounded in L°°(]R d ), uniformly in L e N*. Let 
ip G Ran('P) such that ||V'||L 2 (R d ) = 1) Cl = VLip, and Cl the L7\L-periodic extension of Ql. Then, 
XlH-L,n l Cl G X l , and it holds 



"0 - XillriViCillH 1 



< \\i> - Vl^Whhri) + WCl - Xl^l,n l ClWh 1 (K d ) 

= ||-0 - ^lV'IIhi^) + llXi (Cl - n L!Arj - Ci) llffi(R d ) 

' L + qiCi 



< Ce 

< Ce~ 5L 

cle 



L - H-L,N L (l\\H, 
r-1 

IKilk pV (r 



-8L 



Nl 
L 



IICi||ff''(R< i 



C e' 5L 



C e 



-5L 



— I \W\H-(R d ) 



Nl 
L 



This yields the estimate 

|(l-<^ 
The remaining estimate 



< C e- &L 



C(L 2 (R d ),H 1 («- d )) 



n a T +TZf < c e - 5L , 



Nl 
L 



is a straightforward consequence of ([35 
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6.4 Proof of Theorem IQ1 

Let us first remark that since — — = Gl G N*, Il.m l (Vp er ) = Ii.G t (^per) is a 7^-periodic function. 
Let 4>l,iPl S -Xl be such that </>z, = xlUl and i/>£ = xlvl with ul,vl € Yl^ l . Then, we have 



(Vpcr - Zl,M l (Vpcr)) ^iV'i 



< 



I Vpcr -ZLML(Vper)\ u\ 



1/2 



[ \V peI -l L ,M L (V pe t)\vl 



1/2 



As 



I \V peI -l L ,M L (V p er)\u 2 L = J2 [ \V p er-ll,G L (V P er)\u L (- + R) 

■' Tl Rennr L ^ r 



2 

lt*(r) 



Tl 



< C\\V pci -Ii,g l (V pci )\\ L 2 r (r)||0L||Hi(K d )ll'0L||Hi(R £ i), (36) 



<\\V pcl .-l hGL (V pC r)\\ LU{r) £ IM'+^H 

= C||Vp 0r -Ii !Gi (t/ pcr )|| L 2 cr(r) ||w L ||^i cr(rL) 
< C||Vp 0r -Ii !Gi (^ pcr )|| L 2 cr ( r )||0 L ||^ 1(R£i) , 

we obtain 

(Vpcr - 2i,Mi ( V per ))<f>LlpL 

for a constant C independent of L, with 

II Vpcr -Zl,G L (V pe r)\\L ieA V) < CG^ 

Besides, since W € C°(R d ) n H r - 2 (R d ), 

\\W L -l LML (W L )\\ L 2 {rL) < C 

< C[--) \\W\\ H r-* m , 



L 

Ml 



^erll^-(r) r -> 0. (37) 

per \ / i,— VCXD 



Ml 
L 

Ml 



r-2 



ll^ll^- 2 (r L ) 



r-2 



(38) 



and ||W- Wi||io. (R «i ) < \\W\\ L , 



sup sup 



— > 0. Thus, 

L— j-cjo 

|(ol - <il)(4>l, iI>l)\ 



4>l£X l il> L £X L ll0i||//i(M' i )l|V'L||_f/ 1 (R< 1 ) 



L— >oo 



0. 



Together with the results proved in Section [6~T1 this implies (A2), (A3), (Bl) and (B2) are satisfied 
for 71 = (Xx,,ai,,mi). Assumption (A4) is also satisfied for 71 = {Xl, Q>Li mi,), with 2l(v) 
playing the role of a„(-,-) and m„(y) = m„(-,-) = mz,(v)- To obtain the estimates (|26"T) . ([271 
and (|28jl. it remains to prove that 



Sf. <C 



— j + ||W|U~(R*\r £ _o 



-5£ 



L 

Al 



where 



5£ 



sup sup 

VGRanCP), ||VII I ,2 (s d ) = l*i.eXi 



(a L -5 L )(n^ (Rd V^L) 



Using (13"6")) and (|3T|) . we already have for all ip G Ran("P) such that ||^|| r.^imd^ = 1, 

r-2 



(Vp er — Il,M l (^per)) (b^V) 



< C 



Wr 



||<^||ffi(R<<)- 



(39) 
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Besides, using (l38l) . it holds that 



(w L -i L ,M L (w L )){n^ V)<^ 



<C[ — \ \\4>L\\m(M.«-)- 



(40) 



It also follows from (1351) that 



(w L - w) (n£ (Rd V) H 



< 



T L \T L 



T L \T L 



. Tt H 1 (M d ) , 



< C||W|| i «() R «i\ ri _ 1 ) I e" dJj + 
Reasoning as in the proof of (Al) in Section [Ol and using (|35|) . we can prove that 

-8L ( L 



4>L\\H 1 (WL d )- 



Thus, 



(w L - w) (n^ (Rd V) 



< c 



L 

Ml 



< C 



Nr 



\W\\ 



-SL 



L 

~~N~l 



Finally, using (gDJ) and (jUJ), we obtain 

r-2 



S a L <C 



L 



\W\ 



r i _ 1 ) 



-<5L 



which ends the proof of Theorem 15.21 



(41) 



7 Numerical results 

In this section, we present some numerical results obtained with the software Scilab, illustrating 
the a priori estimates given in Theorem 15. II and 15.21 These results have been obtained with d = 1, 
Vpcr(x) = | sin a: |, W(x) = — 2exp(— |x|) and T = (— 7T,7r]. The particular form of these potentials 
enables us to compute the mass and stiffness matrices analytically (and therefore with no numerical 
integration error). The operator A = — A + Vp Cr + W then possesses a discrete simple eigenvalue 
A ss 1.69 located in the spectral gap [a,j3] of the operator A = —A + yp er where a ~ 1.43 and 
j3 ~ 1.84. The reference values for A and the associated eigenvector (considered in our numerical 
study as the limits L, Nl — > oo) are obtained with L ro f = 40 and N lc f = 1400. 

Figure 1 shows a(H L ,N Tef ) n [1, 2] for L = 6,8, 10, 12, 14, 16, 18 and iV re f = 1400. We can see 
that there is no spectral pollution, as predicted by [8] and Proposition l6.il 

The next series of numerical tests confirms the exponential convergence of the supercell method 
with respect to the size of the supercell. We have compared the eigenvalue closest to A and 
the associated eigenvector obtained for different values of L (L — 6, 8, 10, 12, 14, 16, 18) to the 
reference eigenvalue and eigenvector obtained with L = 40, all these calculations being done with 
A^ef = 1400. Figure 2 shows the relative errors on the eigenvalue, and the square of the L 2 
and H 1 norms of the error on the eigenvector. More precisely, for all L £ N* , we consider the 
eigenvector ul of Hl.n io{ associated with the eigenvalue Al of Hl,n iq[ closest to 1.69, and set 
4>l = Xlul, where %l is the unique C 2 function defined by \h = 1 on [— ttL,ttL], \l = on 
\ [— 7r(L + V~L), n(L + VX)], and xl is a sixth degree polynomial on [— n(L + VL), —ttL] and 

'|Az,-A £ . \\ , I \\4>l-4>l 



on [ttL, n(L 
and log 10 ( - 



-VI)]. Figure 2 shows the decay rate of log^Q 



(^).'oi 1 c( i 



cf "L 2 (B.) 



ef N ff'(l) 



J rof II 1.2(E) 

These numerical results show the exponential decay of the error as 
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Spectrum of the discretized operators for different values of L 
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Figure 1: Spectrum of Hl,n ic( in the range [1,2] for different values of L, with iV rc f = 1400. 
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Figure 2: Decay rates of log 10 ( '^j^"* ' ) (Eigenvalue), log 10 ( nr^TFT^ ) ( ErrOT L2 ) and 



ef M i2(K) 



log 10 ) (En '" , IJ I i [,,r <lil!«')vi): values of I. 



J ref Hfflfffi) 
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Eigenvalue error 




Figure 3: Error on the eigenvalue log 10 {\\l,n l ,m l — Al|) as a function of log 10 (iV). 



a function of L, as well as the doubling of the convergence rate of the eigenvalue with respect to 
the convergence rate of the eigenvector. 

The last series of numerical tests aims at testing the effect of numerical integration. For all 
L £ N*, we denote by \l.n l .m l the eigenvalue of Hl,n l ,Ml which is closest to A, by ul,n Li m l 
an associated normalized eigenvector, and by 4>l,n l ,m l = Xlul.n l ,m l (we choose the sign of 
ul.n l .m l in such a way that \\4>l,n l ,m l — 0L||L 2 (R d ) — 0). In the plots below are drawn the errors 
\^l,n l m l - \\(f>L,N L ,M L - 0L||i2(Rd) and \\<f> L - (t>L,N L ,M L \\H 1 (R*) for tne following values: 

• L = 6,8,10,12,14,16,18, 

• N L = NL where N ~ 2,4, 6, 8, 10, 12, 14, 

• M L = ML where M = 56, 112, 224, 448, 

as well as the results obtained with exact integration (M = oo). 



8 Appendix: Banach-Necas-Babuska's Theorem and Strang's 
lemma 

In this appendix, we recall the Banach-Necas-Babuska theorem and the Strang lemma (see e.g. 

EJUS]). 

Theorem 8.1. (Banach-Necas-Babuska) Let W be a Banach space and V a reflexive Banach 
space. Let a £ C(W x V;R) and f £ V . Then the problem 

find u £ W such that , . 

\/v£V, o(«, «) = /(«), ( ' 

is well-posed if and only if 

-i • r \ a ( W J V )\ 

• da > 0, s.i. ml sup - — n — — — > a; 

ivew veV \\w\\w\\v\\v 
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Figure 4: Error on the eigenvector log 10 (\\c/>l,n l ,m l — 4>L\\L 2 (m d )) as a function of log 10 (iV). 
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• e V, (Vio e W, a(w, u) = 0) (u = 0). 
Moreover, the following a priori estimate holds: 

v/eV, ||«|k<-||/||v. (43) 

a 

Lemma 8.1. (Strang) Le£ ms consider the following approximate problem 

( find u n G VF„ such that .. . 

1 V«„ e a„(n n ,n n ) = /„(«„), 

and fe< ms assume that 

• W„ C W and K C V; 

• 3a n > 0, s.i. inf sup J a "f W "' > a„, and dim(W n ) — dim(V n ); 

w n ew„ VneVn ||u; n ||vH|t>n||v 

• the bilinear form a n is bounded on W n x V n . 
Then, the following error estimate holds: 

\\U-U n \\ w < — 1|/ - fn\\c(V n ) 

+ inf 

w n EW, 



(. . \\<A\c{W,V n ) \ 1 1 II , 1 

H )\\u-w n \\w^ sup 
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